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Abstract 

MicroRNAs (miRs) play a pivotal role in a variety of biological processes including stem cell differentiation and function. 
Human foetal femur derived skeletal stem cells (SSCs) display enhanced proliferation and multipotential capacity indicating 
excellent potential as candidates for tissue engineering applications. This study has examined the expression and role of 
miRs in human foetal femur derived SSC differentiation along chondrogenic and osteogenic lineages. Cells isolated from the 
epiphyseal region of the foetal femur expressed higher levels of genes associated with chondrogenesis while cells from the 
foetal femur diaphyseal region expressed higher levels of genes associated with osteogenic differentiation. In addition to 
the difference in osteogenic and chondrogenic gene expression, epiphyseal and diaphyseal cells displayed distinct miRs 
expression profiles. miR-146a was found to be expressed by human foetal femur diaphyseal cells at a significantly enhanced 
level compared to epiphyseal populations and was predicted to target various components of the TGF-p pathway. 
Examination of miR-146a function in foetal femur cells confirmed regulation of protein translation of SMAD2 and SIV1AD3, 
important TGF-p and activin ligands signal transducers following transient overexpression in epiphyseal cells. The down- 
regulation of SMAD2 and SIVIADB following overexpression of miR-146a resulted in an up-regulation of the osteogenesis 
related gene RUNX2 and down-regulation of the chondrogenesis related gene S0X9. The current findings indicate miR-146a 
plays an important role in skeletogenesis through attenuation of SI\/1AD2 and SMAD3 function and provide further insight 
into the role of miRs in human skeletal stem cell differentiation modulation with implications therein for bone reparation. 
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Introduction 

Skeletogensis is a multistep process consisting of mesenchymal 
cell condensation, proliferation, hypertrophic difTerentiation of 
chondrocytes, and finally, mineralization of extracellular matrix by 
osteoblasts [1-3]. The process of skeletogensis is orchestrated by 
various factors including transcription factors [4], micro environ- 
mental signals and epigenetic cues [5,6]. Defects in the regulators 
of skeletogensis results in skeletal dysplasias, growth failure [2] . A 
clearer understanding of skeletal stem and bone cell formation and 
function is critical to inform bone formation strategies and 
subsequently restore the function of the skeletal system. The cell 
responsible for bone formation, the osteoblast, is derived from a 
multipotential marrow stromal stem cell termed the mesenchymal 
stem cell (undifferentiated multipotent cells of the mesenchyme) 
which has gained wide acceptance, however this term is non- 
specific and the term skeletal stem cell (SSC) will be applied to 
restrict description to stem cells from bone able to generate all 
skeletal tissues. 



MicroRNAs (miRs) are a class of non-protein coding small 
RNA molecules of 21-25 nucleotides in length. Along with the 
RNA-induce-sUencing complex (RICS), they possess the ability to 
regulate protein translation by inhibiting their target mRNAs 
function [7] . There are cumulative evidences to suggest miRs plays 
an important role in many cellular processes including cell cycle 
and stem cell differentiation [8,9]. Various miRs have already 
been identified to play a role in SSC differentiation, a recent 
review by Lian et al have summarized the effects of 42 miRs on 
osteoblast dififerentiation through targeting various cells signaling 
pathways such as Wnt and TGF-P, transcription factors such as 
RUNX2 and Osterix and epigenetic machineries such as histone 
deacetylase 5 (HDAC5) [10]. Data gathered through proteomic 
approach have demonstrated that a single miR can repress the 
production of hundreds of proteins, however, the effect of a single 
miR on protein translation is surprisingly small [11], therefore it 
can be difficult to determine how a single miR is able to provoke a 
detectable functional change. 

Human foetal femur derived SSC have been shown to contain 
stromal antigens positive cells with the potential to differentiate 



PLOS ONE I www.plosone.org 



1 



June 2014 | Volume 9 | Issue 6 | e98063 



MiRNA-146a and Skeletal Stem Cell Differentiation 



down osteogenic, chondrogenic and adipogenic lineages when 
treated with appropriate culture conditions [12]. Furthermore, 
foetal femur cell populations have been shown to possess enhanced 
renewing, differentiation and immunoprivilege potentials, indicat- 
ing their potential as a cell source for tissue engineering 
applications [12,13]. However, cells isolated from the foetal femur 
comprise a heterogeneous population of cells with differing aflBnity 
and capacity for chondrogenic and osteogenic differentiation [14] 
all of which has serv ed to limit their clinical translation. 

A number of growth fac:tors, signaling molecules and transcrip- 
tion factors have been shown to affect skeletal stem cell and 
osteoprogenitor cell activity including members of the Wnt and 
TGF-P families [15-17]. Furthermore, a number of miRs have 
been reported to be involved in the regulation of osteogenesis and 
chondrogenesis through their ability to regulate transcription 
factors [18]. Thus, miR-140 has been identified as a cartilage 
specific miR capable of promoting chondrogenic differentiation by 
increasing the expression of RUNX2, a g(;n(; important in 
chondrocyte hypertrophic differentiation, through down-regulat- 
ing HDAC4 [19,20]. More recently, miR- 138 was reported by 
Eskilden and coworkers to be a negative regulator of osteogenic 
differentiation through inhibition of the expression of Osterix 
(OSX) via targeting focal adhesion kinase (FAK) [21]. Thus an 
understanding of the interactions of specific miRs with signaling 
pathways and growth factors that modulate bone cell fianction 
offers new strategies to manipulate and modulate SSC differen- 
tiation enhancing our understanding of bone physiology and 
function critical in any reparative approach. 

The current study has examined the expression of miRs within 
human foetal femur derived diaphyseal and epiphyseal popula- 
tions. Following identification of select miRs, function was 
examined using transient overexpression analysis for chondrogenic 
and osteogenic differentiation. We demonstrate that cells isolated 
from the epiphyseal regions of the developing foetal femur 
expressed higher levels of chondrogenic related genes while cells 
from the diaphyseal region expressed genes associated with 
osteogenesis. Using RT-qPCR methods, we have confirmed the 
expression of stromal antigens CD63, ALCAM and putative stem 
cell marker Nucleostemin by both epiphyseal and diaphyseal cells. 
MicroRNA microarray analysis confirmed differential miR 
expression in epiphyseal and diaphyseal populations. Furthermore, 
putative target and functional analysis demonstrated that miR- 
146a affected skeletal cell differentiation by down-regulation of 
SMAD2 and SMAD3 protein translation, genes known to be 
involved in the activation of chondrogenesis [22]. The current 
findings indicate an important role of miR- 146a, in conjunction 
with other miRs already described in literature during skeletogen- 
esis and wiU inform our understanding of bone development and 
reparation. 

Method 

Ethics Statement 

Human foetal tissue was obtained with informed and written 
consent from women undergoing termination of pregnancy 
procedure according to guidelines issued by the PoUdnghome 
Report. Ethical approval was obtained from the Southampton & 
South West Hampshire Local Research Ethics Committee (LREC 
296100). 

Isolation and Culture of Foetal Fennur Derived Cells 

Foetal femurs samples at 7-9 weeks post conception were 
isolated from fetuses. Surrounding skeletal muscle and connective 
tissues were removed from the collected foetal femur samples and 



separated into epiphysis (region containing densely packed cells) 
and diaphysis (mid shaft region containing hypertrophic cells) 
sections by micro-dissection. Femurs were plated into a well of a 6- 
well-plate (Costar) overnight in a-MEM (Gibco) containing 1 mg/ 
ml coUagenase B (Roche). The cell suspension was passed through 
a 70 |J,m filter, centrifuged and resuspended in ot-MEM supple- 
mented with 10% foetal calf serum (FCS) and 1% penicillin/ 
streptomycin mix (P/S). Cells were maintained in monolayer 
culture under standard condition until 9Q% cell confluence was 
reached before passaging using fx Trypsin solution (Lonza). All 
monolayer experiments were conducted on passage two cells. 
Foetal age was determined by measuring foetal foot length and 
samples classified as weeks post conception. 

Osteogenic and Chondrogenic Conditions 

Cells were passaged into six-well plates for osteogenic and 
chondrogenic differentiation. Control cultures were refreshed with 
basal medium (a-MEM supplemented 10% FCS) every 48 hours. 
For osteogenic culture, cells were cultured using medium 
supplemented with 10% FCS, 10 nM of dexamethasone and 
100 |a,M of ascorbic acid 2-phosphate. For Chondrogenic culture: 
a-MEM containing 5 ng/ml TGF-PS, 10 nM dexamethasone, 
100 |i.M ascorbate, and 10 |Xl/ml of lOOx ITS liquid media (Sigma 
13146) was used. 

Alcian Blue/Sirius Red Staining 

Whole foetal femurs were fixed with 4% PFA-PBS overnight at 
4°C and processed through graded ethanol. Samples were 
immersed in 50%) chloroform/ethanol and 100"/() chloroform. 
Samples were placed in paraffm wax at 60°C for 30 minutes to 
enable full penetration of paraffin wax. Thereafter, tissues were 
embedded in paraffin wax blocks for sectioning. Sections were cut 
on a Microm 330 at 6 \lm thick and the cut sections were 
transferred onto pre-heated glass slide. Prior to staining, the 
samples were processed through histoclear solution twice to 
remove paraffin wax and rehydrated using reverse graded 
metiianol solutions (100%, 100%, 90% and 50%). To allow 
visualization of the cell nucleus, Weigert's haematoxylin solution 
was prepared and added to rehydrated samples for 10 minutes. 
Excess staining solution was removed with water and acid alcohol. 
To stain for proteoglycans, samples were immersed in 0.5% Alcian 
blue 8GX solution. Samples were placed in 1% molybdopho- 
sphoric acid followed by incubation in 0.1% Sirius red FSB 
solution to stain for collagen. Sections were rinsed thoroughly with 
water and dehydrated in reverse graded methanol back into 
histoclear before mounting in dibutyl phthalate xylene (DPX). 

Image Capture and Analysis 

Sample images were captured using a Zeiss Axiovert 200 
inverted microscope and Zeiss Axiovision software version 4.7. 
Light microscopy images were taken using an Axiocam HR 
camera and fluorescent images were captured using Axiocam MR. 

RNA Extraction and RT-PCR Analysis 

RNA extraction was performed using mirVana RNA Isolation 
System Kit (Life Technologies) according to the manufacturer's 
protocol. Briefly, cultured samples were placed on ice and washed 
twice in PBS. Lysis buffer was used to release RNA and a miR 
homogenizing agent was added, followed by acid phenol- 
chloroform. The resultant mixture was centrifuged to allow phase 
separation. The aqueous phase was removed and added to ethanol 
prior to elution through a spin column. The column was washed 
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Table 1. Forward and reverse primers sequence used for RT-qPCR. 



Gene 


Primer sequences 


Amplicon size 


Human fi- Actin 


F: 5' gge ate etc ace etg aag ta 3' 


82 bp 


R: 5' agg tgt ggt gcc aga ttt tc 3' 


Human Runx-2 


F: 5' gta gat gga eet egg gaa cc 3' 


78 bp 


R: 5' gag gcg gtc aga gaa caa ac 3' 


Human ALP 


F: 5' gga act eet gae cct tga cc 3' 


86 bp 


R: 5' tec tgt tea get cgt act gc 3' 


Human Coll a 1 


F: 5' gag tgc tgt ccc gtc tgc 3' 


52 bp 


R: 5' ttt ctt ggt egg tgg gtg 3' 


Human Type X Collagen 


F: 5' eec act ace caa eae caa ga 3' 


95 bp 


R: 5' gtg gae cag gag tae ctt ge 3' 


Human Osteonectin 


F: 5' gag gaa ace gaa gag gag g 3' 


95 bp 


R: 5' ggg gtg ttg tte tea tee ag 3' 


Human SMAD3 


F: 5' tga ate cct acc act acc aga g 3' 


117 bp 


R: 5' gga tgg aat gge tgt agt eg 3' 


Human SMAD 2 


F: 5' gat cct aae aga act tee gee 3' 


146 bp 


R: 5' cae ttg ttt ete eat ett eae tg 3' 


Human Sox9 


F: 5' ccc tte aae etc cea cac ta 3' 


74 bp 


R: 5' tgg tgg teg gtg tag teg ta 3' 


Human Col2a1 


F: 5' eet ggt ccc cct ggt ctt gg 3' 


58 bp 


R: 5' cat caa ate etc cag cea tc 3' 


Human Nucleostemin 


F: 5' ggg aag ata ace aag cgt gtg 3' 


98 bp 


R: 5' cct cea aga agt tte caa agg 3' 


Human CD63 


F: 5' gcc ctt gga att get ttt gtc g 3' 


87 bp 


R: 5' cat eae etc gta gee act tet g 3' 


Human ALCAM 


F: 5' aeg atg agg cag acg aga taa gt 3' 


96 bp 


R: 5' cag caa gga gga gae caa caa e 3' 



doi:l 0.1 371/journal.pone.0098063.t001 



three times with supplied bufTer solutions and RNA eluted in 
RNase free water. 

For cDNA synthesis, Superscript VILO cDNA Synthesis Kit 
(Life Technologies) was used. RNA was combined with 5X VILO 
reaction mix and 1 |J,1 of lOX Superscript enzyme and incubated 
at 25°C for 10 minutes followed by 42°C for 2 hours. The reaction 
was terminated by incubation at 85°C for 5 minutes. 40 [Xl of 
water was then added to the cDNA sample to give a 1:4 dilution 
and stored at — 20°C or used immediately for quantitative RT- 
PCR analysis. Quantitative RT-PCR was performed using SYBR- 
Green PGR master mix (Life Technologies): 10 |.tl of SYBR-Green 
master mix; 5 |J,1 of up H2O; 2 |Xl of forward and reverse primers 
for the gene of interest (Table 1) and 1 |J,1 of cDNA sample. The 
final mixture (20 |a,l) was added to each well of a 96-weU-plate, 
analyzed using an Applied Biosystem (Life Technologies), 7500 
Real Time PGR system. Data was analyzed using the Applied 
Biosystem 7500 System SDS Software, version 2.0.5. Gt value for 
each sample was normalized to P-Actin, an endogenous house- 
keeping gene, and fold expression levels for each target gene were 
calculated using the delta-delta Gt (Gycle threshold) method. 

MicroRNA Microarray Analysis 

Applied Biosystems TaqMan Low Density Array system was 
used to perform RT-PGR miR arrays, following the manufactur- 
er's protocol. In brief, 400 ng of total RNA was used to generate 
cDNA using MegaplexRT Human Pool A (Life Technologies). 



6 I.1I of cDNA was then combined with 450 |a,l of 2x TaqMan PGR 
Master mix (Life Technologies) and 444 |J,1 of water before loading 
into each micro array card. An Applied Biosystems 7900 
thermocycler was used to run the array. The raw data generated 
were analysed using Sequence Detection System v2.4 Enterprise 
edition (Applied Biosystems). Further data analysis and heat maps 
were created in R (v 2.15.3). Raw expression values were pre- 
processed using the R package qpcrNorm [23]. Heat maps were 
generated from quantHe normalized, z-transformed Gt values, 
ordered by ascending ddGt of epiphysis versus diaphysis. A subset 
of all expressed miRss was selected for a heat map based on 
statistical evidence for differential expression according to a 
permissive t-test (o(=0.1). 

MicroRNA Expression Analysis 

MicroRNA expression was examined using TaqMan Micro- 
RNA Assays (Life Technology). Each assay comprised two 
primers, one for cDNA synthesis and one for TaqMan RT- 
PGR. cDNA specific to each assay-specific miR was generated 
using TaqMan MicroRNA Reverse Transcription Kit (Life 
Technologies) from total RNA using a modified manufacturer's 
protocol. In brief, a reaction mixture containing 50 mM dNTPs, 
25 units of Mutiscribe RT enzyme, 0.75 |il lOx Buffer, 1.88 units 
of RNase inhibitor, 3.58 |J,1 water, 1.5 |J,1 of primer and 10 ng of 
total RNA was prepared. The reaction mixture was mixed by 
pipetting and incubated at 16°G for 30 minutes followed by 42°G 
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Figure 1. Morphology of foetal femur derived cells in whole femur and post monolayer culture with quantitative expression of 
stromal antigens and putative stem cell marker by HBMSC, foetal femur epiphyseal and diaphyseal cells. Foetal femur comprised a 
proteoglycan aniage (Blue) with a bone collar marked by deposition of collagen (Red) (A). Epiphyseal region contained proliferating chondrocytes (B) 
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while diaphyseal section, comprised hypertrophic chondrocytes (C). Epiphyseal cells (D) and diaphyseal cells (E) adopt similar morphology to HBMSC 
(F) upon monolayer culture. Expression of CD63 (G), ALCAM/CD166 (H) and Nucleostemin (I) by foetal femur epiphyseal cells and diaphyseal cells was 
confirmed by RT-qPCR and expression levels compared to human bone marrow stromal cells (set to have expression of one). Results expressed as 
mean ± SD and n = 3. ***P<0.001 and *P<0.1 calculated using ANOVA. Scale bar= 100 nm. 
doi:l 0.1 371 /journal.pone.0098063.g001 



for 30 minutes. The reaction was terminated following incubation 
for 5 minutes at 85°C. The cDNA sample was stored at — 20°C if 
not used immediately for qPCR analysis. 

For the qPCR reaction, TaqMan Universal PGR Master Mix 
with No AmpErase UNG (Life Tc-clniologics) was utilised. A 
reaction mixture containing 5 \il TaqMan master mix, 3.335 |Xl of 
up H2O, 0.5 Hi of qPCR primer and 0.7 Hi of cDNA was 
prepared. The mixture was transferred into the well of a 96-weU 
plate and loaded into the Applied Biosystem (Life Technology, 
USA), 7500 Real Time PGR system for assay. All reactions were 
performed in duplicate and included a negative control lacking 
cDNA. The data was analyzed using the Applied Biosystem 7500 
System SDS Software, version 2.0.5 program. Gt value for each 
sample was normalized to MammU6, an endogenous control for 
miR and fold cxprc-ssion levels for each target gene was calculated 
using the delta-delta Gt (Gycle threshold) method. 

MicroRNA Transient Overexpression 

To assess the effect of miR- 146a on cultured cells, transient 
over-expression of miR- 146a was achieved using a miR- 146a 
mimic and DharmaFECT reagent (Dharmacon). Briefly, epiph- 



yseal cells were expanded in culture until 80% confluent. Gells 
were trypsinised and plated to each well of a 24-weU plate in ot- 
MEM supplemented widi 10% FGS and 1% P/S and cultured for 
24 hours. On transfection, cells were washed with PBS and 
transfection media containing DharamFEGT reagent and 100 r|m 
of miR- 146a mimic was added to cultured cells. For the control 
group, [:ells were transfected with a scrambled miR mimic 
(Dharmacon) under the same condition. RNA was then extracted 
for RT-PCR analysis 48 hours and protein was obtained after 72 
hours post transfection for western blot analysis. 

Western Blot Analysis 

Gells were lysed with RIPA buffer (750 mM NaGl, 5% IgePal 
GA-630, 2.5% DOG, 0.5% SDS, 250 mM Tris (pH 8.0)) 
supplemented with a protease inhibitor cocktail (1:100; Sigma). 
Gell lysates were centrifugation at 13000 rpm for 20 minutes at 
4°G, and the supernatant collected. Protein concentration was 
determined using Pierce BCA protein Assay Kit (Pierce). 20 |Xg of 
each samples combined with DTT (GeU Signalling) and SDA 
sample buffer (BioLabs) were analysed by SDS gel electrophoresis 
using Any kD precast gels (Mini-PROTEAN TGX; Bio-rad) and 
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Figure 2. Relative expression (RT-qPCR) of osteogenic and chondrogenic marker genes in epiphyseal and diaphyseal cell 
population following monolayer culture. Cells extracted from the diaphyseal region expressed higher levels of osteogenic marker genes; 
RUNX2, ALP, Type I collagen and Osteonectin (A-D) while epiphyseal cells expressed higher level of genes associated with chondrogenesis; S0X9 and 
Type II collagen (E-F). Relative expression was normalized to P-Actin with epiphyseal cell populations set to have an expression of one. Data are 
represented as an average of three independent patient samples and error bars represent standard deviation. **P<0.01; ***P<0.001. 
doi:l 0.1 371/journal.pone.0098063.g002 
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Figure 3. RT-PCR results showing tKie effect of osteogenic media on osteogenic and cKiondrogenic gene expression in epiphyseal 
and diaphyseal cell populations. Epiphyseal cells, in osteogenic media, expressed higher levels of RUNX2, ALP and Type I collagen (A-C) whilst in 
diaphyseal cell populations, only an increase in ALP expression was recorded (E-G). In epiphyseal and diaphyseal cells, the expression of the 
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chondrogenic transcription factor, 50X9, was reduced in osteogenic media (D and H). Relative expression was normalized to P-Actin with basal 
condition set to an expression of one. Data represent an average of three independent patient samples and error bars represent standard deviation. 
**P<0.005; ***P<0.001 calculated using Mann-Whitney test. 
doi:1 0.1 371/journal.pone.0098063.g003 



transferred onto a polyvinylidene fluoride (PVDF) immobilon-FL 
transfer membrane (Millipore). Immunoblots were probed with 
SMAD3 and SMAD2 antibody (1:1000; CeUsignalling), then with 
anti-P-Actin-Peroxidase (1:10000; Sigma). IRDye 800CW goat 
anti-rabbit (LI-COR Biosciences) and IRDye 800CW rabbit anti- 
HRP (LI-COR Biosciences) were used as secondary fluorescent 
antibodies respectively. The anti-P-actin-Peroxidase was used as 
loading control to which protein expression across the membrane 
was normahzed. The blots were scanned individually on an 
Odyssey Infrared Imaging System (LI-COR Biosciences), and 
densitometry analyses were run calculating the signal intensity 
ratio between samples and loading control. Image Studio software 
(LI-COR, Biosciences) was employed to obtain the signal intensity 
values which were transformed into percentages of intensity for the 
statistical analysis. 

Statistics 

Statistical analysis was carried out using Mann-Whitney test or 
One-way Analysis of Variance (ANOVA) with Tukey-Kramer 
multiple comparison post-test using statistics software package 
Prism version 6 (GraphPad software). Values were expressed as 
mean ± standard deviation. AH experiments were performed 
using at least three separate; foetal samples unless otherwise stated. 
Values for p£0.05 were considered statistically significance. 

Results 

Isolation of Osteogenic Diaphyseal and Chondrogenic 
Epiphyseal Foetal Femur Cells Population 

Foetal femurs were observed to comprise primarily of a cartilage 
anlage with a developing bone collar marked by collagen 
deposition, as evidenced by Alcian blue and Sirus red staining 
(Figure lA). Within the epiphyseal regions, cells were observed to 
be smaller and densely packed (Figure IB), while cells from the 
diaphyseal region displayed a hypertrophic phenotype (Figure IC). 
Following monolayer culture, both epiphyseal and diaphyseal 
derived cell populations formed distinct cell colonies displaying a 
fibroblastic morphology similar to the human bone marrow 
stromal cells (Figure ID-F). RT-qPCR confirmed expression of 
stromal antigens CD63 and ALCAM (Figure IG— H) in epiphyseal 
and diaphyseal cells although at a lower level in comparison to 
human bone marrow stromal cells (HBMSC) (Figure IG). 
Epiphyseal cells expressed ALCAM at a similar level to HBMSC 
while diaphyseal cells appeared to express lower levels of ALCAM 
(Figure IH). Nucleostemin, a putative stem cell marker, was found 
to be expressed by both epiphyseal and diaphyseal ceUs (Figure II). 
The expression of CD63, ALCAM and Nucleostemin by 
epiphyseal and diaphyseal cells in monolayer culture, suggest 
these cell populations retain mesenchymal progenitor cell-like 
characteristics. We have previously published on the potential of 
foetal femur derived cells to differentiated along the stromal 
osteogenic, chondrogenic and adipogenic lineages [12] and this 
was not re-examined here. 

Epiphyseal and Diaphyseal Cells Display Distinct 
Osteogenic and Chondrogenic Differentiation Potential 

Epiphyseal and diaphyseal cells were isolated from three 
unrelated patient samples. Following 7 days in monolayer culture 
under basal condition, RNA was extracted from each individual 



sample for gene expression analysis. Using RT-qPCR, diaphyseal 
cells were shown to express higher levels of genes associated with 
the osteoblast phenotype, namely RUNX2, ALP, Type I Collagen 
and Osteonectin (Figure 2A-D). Epiphyseal cells expressed 
increased level of genes associated with the chondrocyte pheno- 
type, namely SOX9 and type II collagen (Figure 2E-F). Following 
culture in osteogenic conditions, epiphyseal and diaphyseal cells 
expressed higher levels of ALP mRNA (Figure 3B and F), while, 
RUNX2 and Type I Collagen expression were only observed to be 
increased in the epiphyseal cell population (Figure 3A and C), 
suggesting (enhanced osteogenic modulation of epiphyseal cell 
populations by osteogenic media. In epiphyseal and diaphyseal 
cells, SOX9, a gene associated with chondrocyte differentiation, 
was reduced following osteogenic media supplementation 
(Figure 3D and H). 

Epiphyseal and Diaphyseal Cells Express Distinct 
MicroRNAs 

RT-qPCR miR array was used to examine the miR expression 
profile in the chondrogenic epiphyseal cells and osteogenic 

diaphyseal cells in three unrelated foetal femur samples. The 
miR array expression analysis indicated 1,5.^ out of the 377 
detectable miRs were' cxprc'sscxl in the samples. Of tlu^sc", 67 miRs 
were found have a difference in expression of greater than 1.5 fold 
between epiphyseal and diaphyseal cells with 1 2 miRs expressed at 
an elevated level in epiphyseal cells and 55 miRs in the diaphyseal 
cell population. Only 7 miRs were identified with a statistically 
significant difference in expression between epiphyseal and 
diaphyseal cells (Figure 4A), namely: miR-146b-5p, miR-301b 
and miR- 138 (higher expression in epiphyseal cells) and miR- 143, 
miR-145, miR-146a and miR-34a (increased expression in 
diaphyseal cells). A heat map was employed to demonstrate the 
difference in selected miRs expression between epiphyseal and 
diaphyseal cells (Figure 4B). miR- 146a was found to be expressed 
at a higher level in the diaphyseal cell population and was selected 
for revalidation using an individual TaqMan RT-qPCR assay. 
miR- 138 and miR- 140, have previously been reported to have 
anti-osteogenic and pro-chondrogenic properties respectively 
[21,24] and were found to display a higher expression in 
epiphyseal cell populations and thus selected for revalidation 
using individual TaqMan RT-qPCR assay to assess the consis- 
tency of data of the current study compared to current literature. 
In addition, the effects of osteogenic and chondrogenic media on 
the expression of miR- 140, miR- 138 and miR- 146a were 
examined. Individual TaqMan assays confirmed the expression 
of the cartilage specific miR- 140 and the anti-osteogenic miR- 138 
was higher in the epiphyseal cell populations (Figure 5A and B). 
miR- 146a displayed a 50-fold increase in expression in diaphyseal 
cells relative to epiphyseal cells (Figure 5C). Culture in osteogenic 
media failed to modulate the expression of miR- 138, miR- 140 and 
miR- 146a (Figure 5D and E). Critically, the expression of miR- 
146a was markedly reduced in the presence of chondrogenic 
media (Figure 5F). 

MicroRNA-146a Target Analysis 

A number of the miRs identified by the microarray analysis 
have already been described to be involved in the differentiation of 
SSC, namely; miR-146b-5p [25,26], miR-138 [7,21], miR-143 
[8,9,27] and miR-145 [10,28,29]. However, the effects of miR- 
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Figure 4. MicroRNA expression profile of epipliyseal and diapKiyseal cells. CT values of each mlR was normalized to MammU6 (dCT) and 
plotted as an XY scattered chart displaying the correlation of mlR expression by epiphyseal and diaphyseal cells. An unpaired t-test revealed 7 mlRs 
with significant differences in expression: miR-146a-5p, miR-301b and miR-138 displayed higher expression in epiphyseal cells while miR-143, mlR- 
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146a, miR-34a and miR-145 displayed higher expression levels in diaphyseal cells (A). A heat map of normalized CT values was generated to show the 
difference in miR expression between epiphyseal and diaphyseal cell populations (B). The Heat map represents z-transformed expression, blue 
represents higher level of expression and red represents lower level of expression. 
doi:10.1371/journal.pone.0098063.g004 



146a on bone cell function are unknown. Furthermore, through 
putative target analysis using Targetscan [30] , miR- 146a appeared 
to target various key regulators the TGF-P ligand specific pathway; 
SMAD2, SMAD3 and SMAD4. SMAD2 and SMAD3 were 
chosen for validation as they are known intracellular transducers 
of TGF-P and activin ligands and are thought to play a critical role 
in chondrogenic differentiation [11,22] whilst SMAD4 was 
omitted in the analysis as the effect of miR- 146a on SMAD4 
function have already been described [12,31]. The putative 
binding site of miR-146a to SMAD2 and SMAD3 3'UTR is 
shown in figure 6A and 6B respectively. To determine whether a 
correlation existed between miR-146a expression and SMAD2/ 
SMAD3 expression in epiphyseal and diaphyseal cell populations, 
RT-qPCR expression was undertaken. The increased expression 



of miR- 146a observed in diaphyseal cells was associated with 
decreased expression of SMAD3 mRNA level however; no 
correlation to the expression of SMAD2 was observed 
(Figure 6C-E). Furthermore, enhanced expression of miR- 146a 
was associated with an increased expression of RUNX2 and ALP 
(6F-G). Critically, decreased expression of SOX9 was observed 
when miR- 146a was highly expressed (Figure 6H). 

MicroRNA-146a Regulates Skeletal Stem Cell 
Differentiation by Down-regulating SMAD2 and SMAD3 
Protein Translation 

To validate the effect of miR- 146a on SMAD2 and SMAD3 
protein translation, miR- 146a was transiendy overexpressed in 
epiphyseal cells from three independent patient samples using a 
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Figure 5. Relative expression (RT-qPCR) of miR-140, miR-138 and miR-146a in epipKiyseal and diaphyseal cell populations (A-C) and 
the effect of differentiation media on expression in diaphyseal cells (D-F). miR-140 was expressed at a higher level in epiphyseal cells (A). 
miR-138, known to have anti-osteogenic effects was found to have a lower level of expression in diaphyseal cells (B). miR-146a, identified in our 
microarray, was validated and shown to display a markedly increased expression level in the diaphyseal cell populations compared to epiphyseal cells 
(C). The use of osteogenic media did not affect the expression levels of miR-140, miR-138 and miR-146a compared to culture in basal conditions (D- 
F). Chondrogenic culture conditions resulted in increased expression of miR-1 40 while miR-1 46a expression was suppressed (D-F). Relative expression 
was normalized to MAMIVI-U6 with epiphyseal cells (A-C) and basal condition (D-F) set to an expression of one. Data represent an average of three 
independent patient samples and error bars represent standard deviation. *P<0.05; **P<0.005; ***P<0.001 calculated using Mann-Whitney (A-C) 
and ANOVA test (D-F). 
doi:1 0.1 371/journal.pone.0098063.g005 
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Figure 6. Predicted binding site for miR-146a In SMAD2 and SMAD3 mRNA 3'UTR and the relationsKiip between expression of miR- 
146a, SMAD2, SIVIAD3 and Type X Collagen expression in epiphyseal and diaphyseal cell populations. SMAD2 and SMAD3 mRNA 
3'UTR contains a binding site for miR predicted by TargetScan (A-B). Increased expression of miR-146a in diaphyseal cell populations was correlated 
with reduced expression in SMAD3 mRNA levels but not SMAD2 mRNA levels (D-E). Increased expression of miR-146a in diaphyseal cells was coupled 
with an increase in expression of RUNX2 and ALP (F-G). Reduced expression of S0X9 was observed when miR-146a was highly expressed. Relative 
expression was normalized to MAMM-U6 (C) and |3-Actin (D-H). Data represents an average of three independent patient samples and error bars 
represent standard deviation. *P<0.05, **P<0.01; ***P<0.001 calculated using Mann-Whitney test. 
doi:10.1371/journal.pone.0098063.g006 
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Figure 7. Effect of miR-146a overexpression on SMAD2/SIVIAD3 mRNA expression and protein translation and in osteogenic related 
gene expression. Overexpression of miR-146a resulted in a reduction of SMAD3 mRNA expression while SMAD2 mRNA levels remained unchanged 
(A-B). Protein level of both SMAD2 and SMAD3 were significantly reduced by overexpression of miR-146a as evidenced by western blot analysis (C) 
with densitometry band intensity quantification (D-E). miR-146a overexpression resulted in increased RUNX2 expression and a reduction in 50X9 
expression (F-G). Relative expression analysis by RT-qPCR was normalized to p-Actin. Densitometry data presented as percentage of band intensity 
compared to control group (100%). Data are represented as an average of three independent patient samples and error bar represents ± SD. *P< 
0.05, **P<0.01; ***P<0.001 calculated using Mann-Whitney test. 
doi:l 0.1 371/journal.pone.0098063.g007 



miR-146a mimic. RNA was extracted 48 hours post transfection 
and protein was harvested after 72 hours. These time points were 
chosen after testing the transfection agents used in a series of 
optimization experiments (data not shown). 48 hours post miR- 
146a over-expression; down-regulation of SMAD3 mRNA was 
observed however, SMAD2 mRNA levels remained unchanged 
(Figure 7A-B). Western Blot analysis showed a reduction in both 
SMAD2 and SMAD3 protein translation in the presence of miR- 
146a overexpression. Quantitative analysis using western blots 
demonstrated miR-146a significantiy reduced SMAD3 protein 
expression over SMAD2 expression (65% and 35% respectively) 
(Figure 7D-E). The down-regulation in SMAD2 and SMAD3 was 
coupled with a down-regulation of SOX9 expression and up- 
regulation of RUNX2 mRNA expression suggesting miR-146a has 
a positive effect on osteogenesis (Figure 7F-G). 

MicroRNA-146aArGF-p Feedback Mechanism Regulates 
Chondrocyte Hypertrophic Differentiation 

Human foetal cells extracted from the epiphyseal layers undergo 
hypertrophic differentiation in the presence of chondrogenic 
media containing TGF-P3 (Figure 8 A-B). In vitro stimulation of 
epiphyseal cells with chondrogenic media containing TGF- P3 
resulted in down-regulation of miR-146a (figure 8C) and a 
substantial up-regulation of Type X collagen mRNA expression 
(Figure 8D). To further validate the effect of miR-146a as a 
negative regulator of TGF-P ligand dependent signaling pathway 
through down-regulating protein translation of SMAD2 and 
SMAD3, the effect of TGF-P3 on Type X Collagen in the 
presence of miR- 146a overexpression in basal and chondrogenic 
conditions was compared. Under basal condition, overexpression 
of miR- 146a had no effect on the expression of Type X Collagen 
(Figure 8E); however, in the presence of miR- 146a overexpression, 
up-regulation of Type X Collagen by TGF-P3 was reduced by 
60%. These results suggesting the effects of miR-146a on the 
TGF-p pathway are dependent on the presence TGF-PS ligands 
(Figure 8F). 

Discussion 

MicroRNAs have recentiy been demonstrated as important 
regulators of a variety of biological processes [12,13,32] including 
cell cycle, oncogenesis and stem cell differentiation [9,14,33,34]. 
Using TaqMan RT-rjPCR MicroRNA array, epiphyseal and 
diaphyseal cells were found to express a different repertoire of 
miRs. Microarray data are known to be relatively noisy [15- 
17,35] and together with the accepted variation inherent to 
primary cells samples, any statistical analysis of microarray data 
can prove challenging. miR- 146a was identified to be expressed at 
a higher level in diaphyseal cells compared to epiphyseal cells, 
suggesting a role of miR- 146a in osteogenic differentiation. Using 
Targetscan release 6, potential mRNA targets of miR- 146a were 
identified [30]. Various components of the TGF-P pathway; 
namely SMAD2, SMAD3, SMAD4, TGFB-induced factor 
homeobox 1 (TGIFl), BMP and activin membrane-bound 
inhibitor homolog (BAMBI) and activin A receptor type IC/I/ 
Type Il-like 1 (ACVRIC/ACVRIB/ACVRLI) were highlighted 



as potential targets of miR- 146a, advocating the function of this 
miR in diaphyseal cells maybe mediated via attenuation of the 
TGF-P pathway [30]. miR- 146 has previously been identified to 
modulate myofibroblast trans-dififerentiation during TGF-pi 
induction by targeting SMAD4 [18,36] and miR-146 may also 
be an important regulator during the inflammatory state of 
osteoarthritis, as IL- 1 p induced production of TNF-a, a pro- 
inflammatory cytokine known to play a role in osteoarthritis, was 
significantiy reduced by miR-146 overexpression [19,20,37]. 
Furthermore, overexpression of miR-146a has been shown to 
protect the human bronchial epithelial from apoptosis and to 
promote cell proliferation through up-regulation of Bcl-XL and 
STAT3 phosphorylation [21,38]. 

miR- 146a transient overexpression in epiphyseal cells (low level 
expression of miR- 146a compared to diaphyseal cells under 
normal culture conditions) for 48 hours resulted in a significant 
down-regulation of SMAD3 at the mRNA level with reduction of 
SMAD2 and SMAD3 protein level at 72 hours, evidenced by 
western bolt analysis. Concomitant with the reduced SMAD2 and 
SMAD3 levels observed, reduced expression of the chondrogenesis 
related gene SOX9 and up-regulation of the osteogenesis related 
gene RUNX2 was observed. These data suggest miR- 146a is a 
negative regulator of chondrogenesis through down-regulation of 
SMAD2/SMAD3 and may indirectly promote osteogenic differ- 
entiation. Furthermore, as over-expression of miR- 146a resulted 
in a significant negative effect on chondrogenesis in culture over a 
48-72-hour period, the effect of miR-146a under normal 
physiological conditions would appear to be significant. To 
optimize our transfection protocol, miR mimics against GAPDH 
were used and shown to reduce GAPDH expression by over 80% 
after 48 hours post transfection. However, miR inhibitors were 
found to have display a much lower efficacy. Coupled with the 
high expression of miR- 146a in diaphyseal cells, it proved difficult 
to reproducibly demonstrate the effects of miR- 146a inhibitor on 
human foetal diaphyseal cells. 

In the current study, we observed that cells extracted from the 
epiphyseal and diaphyseal regions of a developing human foetal 
femur retain their progenitor cell characteristics, evidenced by the 
expression of CD63, ALCAM and Nucleostemin mRNA, but 
exhibit different and distinct phenotypes and display discrete 
affinities in differentiation along the osteogenic and chondrogenic 
lineages. Foetal epiphyseal cells expressed higher levels of miRs 
inc:luding miR- 140 and miR- 138 reported to promote chondro- 
genesis [20-22] while diaphyseal cells expressed miRs including 
miR-210 and miR-93 previously reported to promote osteogenesis 
[3,23,39]. These data suggests the chondrogenic epiphyseal cell 
population expressed increased levels of miRs associated with 
chondrogenesis while the osteogenic diaphyseal cells expressed 
miRs associated with osteogenesis. 

TGF-P signaling is important for skeletogenesis. It is generally 
accepted that the bone morphogenic proteins (BMPs) and their 
receptors induce early cartilage formation and stimulate mesen- 
chymal cells to differentiate into osteoblasts whilst TGF-P ligands 
and their receptors regulate chondrocyte proliferation and 
differentiation [12,40]. In the current study, the effect of TGF- 
P3 stimulation in monolayer culture model were consistent with 
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Figure 8. The effect of chondrogenic media containing 4 ng/ml of TGF-ps on morphology of cells in monolayer culture and the 
expression of miR-146a and Type X Collagen mRNA. TGF-p3 addition (chondrogenic media) culminated in cells with a distinct hypertrophic 
phenotype (A-B). Expression of miR-146a was reduced following culture under chondrogenic conditions (C) and a significant up-regulation of Type X 
collagen was observed (D). Under basal condition, miR-146a expression did not affect the expression of Type X Collagen (E). Under chondrogenic 
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condition, miR-146a over-expression reduced the stimulatory effect of chondrogenic media on type X collagen (F). 

bars represent 100 |im. 

doi:1 0.1 371 /journal.pone.0098063.g008 



= P<0.05, *»* = P<0.001. Scale 



current literature with TGF-PS stimulated cells differentiating to 
give hypertrophic chondrocytes coupled with a significant up- 
regulation of Type X collagen niRNA expression [21,24,40]. 
TGF-Ps signals are known to be transduced to the nuclei by the 
intracellular mediators, SMADs [41]. To date, eight different 
SMAD proteins have been identified, classified into three 
categories based on their functions; the receptor-activated SMADs 
(SMAD- 1, 2, 3, 5 and 8), common mediator SMAD (SMAD4) 
and the inhibitory SMADs (SMAD6 and 7) [41]. The receptor- 
activated SMADs are firrther divided into two groups based on 
their attachment to ligand specific receptors; SMAD2 and 3 
transduce signal by TGF-f5s and activin ligands while SMADl, 5 
and 8 respond to BMPs stimulation [42]. As SMAD2 and SMAD3 
are known to transduce TGF-fi ligand signals [43], the effect of 
miR-146a on the TGF-(5 pathway was examined through analysis 
of the effects of miR-146a transient overexpression on cells 



stimulated by TGF-fi3. In the absence of miR-146a overexpres- 
sion, TGF-fi3 was observed to up-regulate the expression of Type 
X coUagen 2000-fold. However, following transfection with miR- 
146a mimic, TGF-P3 induced up-regulation of collagen X was 
reduced by over 60%. These results suggest miR-146a attenuates 
the TGF-pS ligand signal, and possibly activin signals, through 
down-regulation of SMAD2 and SMADS protein. Interestingly, 
when cells were stimulated with TGF-P3, a reduction of miR-146a 
expression was observed. The current data indicate the presence of 
a negative feedback mechanism between TGF-P3 stimulation and 
miR-146a expression advocating the presence of an auto- 
regulatory mechanism (Figure 9). Thus miR-146a may reduce 
TGF-P signaling and in turn TGF-PS stimulation may suppress 
miR-146a expression. A similar auto-regulatory feedback mech- 
anism has also been described for the regulation of miR-93 and 
downstream Osterix expression during osteoblast mineralization 
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Figure 9. Modulation of skeletal cell differentiation by miR-146a. miR-146a down-regulated SMAD2 and SMADS protein levels, resulting in 
an attenuation of TGF-p signaling following TGF-|33 stimulation. Stimulation by TGF-|33 down-regulated the expression of miR-146a indicating a 
double negative feedback loop and thus a potential auto-regulatory mechanism. Overexpression of miR-146a under basal condition revealed a 
modest, but positive effect on osteogenesis related gene expression and is likely part of a mlR network involved in promoting osteogenesis. The 
inhibitory effect of miR-146a on SIV1AD4 was previously reported by Zhong et al [31] and was not reevaluated in this study. 
doi:1 0.1 371 /journal.pone.0098063.g009 
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[3] as well as for miR-140 on TGF-|3 signaling through SMAD3 
suppression [44]. 

Various reports have already demonstrated the importance of 
miRs during SSC differentiation; miR-140, a cartilage specific 
miR [19], was reported as a positive effector of chondrogenesis 
through PDGF signaling in zebrafish [45]. In addition, miR-140 
has been linked to the regulation of SMAD3 dependent TGFP 
pathway through down regulation of SMAD3 protein levels and 
thus to play a role in chondrocyte development [44]. miR- 138 has 
been shown to inhibit osteogenic differentiation in telomerase 
immortalized bone marrow derived hMSC through down- 
regulation of FAK and subsequentiy down regulation of the 
FAK downstream targets RUNX2 and Osterix [21]. These studies 
demonstrate various miRs may work in concert to regulate the 
complex mechanisms underlying SSC differentiation. Together 
with the distinct expression patterns observed in epiphyseal and 
diaphyseal cell populations, the data in this study suggests, by 
extracting cells from distinct regions of the human foetal femur, a 
more homogenous skeletal stem cell population can be isolated, 
providing an innovative approach for identification of novo miRs 
involved in skeletal stem cell differentiation and skeletogenesis. 

References 

1. Karsenty G (1998) Genetics of skeletogenesis. Dev Genet 22: 301-313. 
doi:10.1002/(SICI)1520-6408(1998)22: 4<301::AID-DVG1>3.0.CO;2-A. 

2. Michigami T (2014) Current Understanding on the Molecular Basis of 
Chondrogenesis. Clin Pediatt Endocrinol 23: 1-8. doi:10.1292/cpc.23.1. 

3. Yang L, Cheng P, Chen C, He H-B, Xie G-Q, et al. (2012) miR-93/Sp7 
function loop mediates ostcohlast mineralization. J Bone Miner Res 27: 1598— 
1606. doi:10.1002/jhmr.l621. 

4. Pratap J, Calindo M, Zaidi SK, Vradii D, Bhat BM, ct al. (2003) Cell growth 
regulatory role of Runx2 during prolilerative expansion of preosteoblasts. 
Cancer Res 63: 5357-5362. 

5. Zimmermann P, Boeuf S, Dickhut A, Boehmer S, Olek S, et al. (2008) 
Correlation of COLlOAl induction during chondrogenesis of mesenchymal 
stem cells with demethylation of two CpG sites in the COLlOAl promoter. 
Arthritis Rheum 58: 2743-2753. doi:10.1002/art.23736. 

6. Suomi S, Taipaleenmaki H, Seppanen A, Ripatti T, Vaananen K, et al. (2008) 
MicroRN.As regulate osteogenesis and chondrogenesis of mouse bone marrow 
stromal cells. Cene Rcgul Syst Bio 2: 177-191. 

7. Pillai RS i200;i) .MicroRNA function: Mukiple mechanisms for a tiny RNA? 
RNAll:17,i3 1761. doi:10.1261/rna.2248605. 

8. Xu N, Papagiannakopoulos T, Pan G, Thomson JA, Kosik KS (2009) 
MicroRNA-145 Regulates OCT4, SOX2, and KLF4 and Represses Pluripo- 
tency in Human Embryonic Stem Cells. CeU 137: 647-658. doi:10.1016/ 
j.ceU.2009.02.038. 

9. Ivanovska I, Cleary MA (2008) Combinatorial microRNAs: working together to 
make a difference. Cell Cycle 7: 3137-3142. 

10. L.an JB, Stein GS, van Wijnen AJ, Stein JL, Hassan .\1Q, et al. (2012) 
MicroRNA control of bone formation and homeostasis. Nat Rev Endocrinol 8: 
212-227. doi:10.1038/nrcndo.201 1.234. 

11. Back D, ViUcn J, Shin C, Camargo FD, Gygi SP, et al. (2008) The impact of 
microRNAs on protein output. Nature 455: 64-71. doi:10.1038/naturc07242. 

12. Mirmalek-Sani S-H, Tare RS, Morgan SM, Roach HI, WUson DI, et al. (2006) 
Characterization and Multipotentiality of Human Fetal Femur-Derived Cells: 
Implications for Skeletal Tissue Regeneration. Stem Cells 24: 1042-1053. 
doi:10.1634/stemcclls.2005-0368. 

13. Tuch B (1988) Immunogenicity of fetal tissue. Hastings Cent Rep 18: 44. 

14. Developmental plasticity of htmian foetal femur-derived cells in pellet culture: 
self assembly of an osteoid shell around a cartilaginous core. (2011) 
Developmental plasticity of human foetal femur-derived cells in pellet culture: 
self assembly of an osteoid shell around a cartilaginous core. Eur Cell Mater 2 1 : 
558-567. 

15. Back WY, Kim JE (201 1) Transcriptional regulation of bone formation. Front 
Biosci (Schol Ed) 3: 126-135. 

16. Zhou S (201 1) TGF-P regulates P-catenin signaling and osteoblast differentiation 
in human mesenchymal stem cells. J Cell Biochem 112: 1651—1660. 
doi:10.1002/jcb.23079. 

17. DAY T, Guo X, GARRETTBEAL L, Yang Y (2005) Wnt/|3-Catenin Signaling 
in Mesenchymal Progenitors Controls Osteoblast and Chondrocyte Diiferenti- 
ation during Vertebrate Skeletogenesis. Developmental Cell 8: 739-750. 
doi:10.1016/j.devcel.2005.03.016. 

18. Dong S, Yang B, Guo H, Kang F (2012) MicroRN.\s regulate osteogenesis and 
chondrogenesis. Biochem Biophys Res Commun 418: 587-591. doi:10.1016/ 
j.bbrc.2012.01.075. 



Conclusion 

This report has further characterized foetal derived SSC found 
in foetal femur by assessing differentiation related gene expression 
and miR expression profile in epiphyseal and diaphyseal cells. 
Using miR ana)', miR assay and transient over-expression 
analysis, we have identified miR- 146a as an important regulator 
of TGF-fi signaling during chondrocyte development and by 
extension, foetal skeletogenesis, through regulation of SMAD2 and 
SMAD3 protein translation. Further functional analysis to define 
the precise role of miR- 146a and other miR targets identified in 
the epiphyseal and diaphyseal cell populations could provide a 
novel strategy to manipulate SCC differentiation in regenerative 
medicine applications. 

Author Contributions 

Conceived and designed the experiments: KSCC ROCO TS. Performed 
the experiments: KSCC NS PSS. Analyzed the data: KSCC NS PSS TS 
ROCO. Contributed reagents/materials/analysis tools: DIW TS ROCO. 
Wrote the paper: KSCC ROCO TS DIW. 



19. Swingler TE, Wheeler G, Carmont V, ElUott HR, Barter MJ, et al. (2012) The 
expression and function of microRNAs in chondrogenesis and osteoarthritis. 
Arthritis Rheum 64: 1909-1919. doi:10.1002/art.34314. 

20. Tuddenham L, Wheeler G, Ntounia-Fousara S, Waters J, Hajihosseini MK, et 
al. (2006) The cartilage specific microRNA-140 targets histone deacetylase 4 in 
mouse cells. FEBS Lett 580: 4214-4217. doi:10.1o'l6/j.lebslct.2006.06.080. 

21. Eskildsen T, Taipaleenmaki H, StcnvangJ, Abdallah BM, Ditzel N, et al. (201 1) 
MicroRNA- 138 regulates osteogenic differentiation of human stromal (mesen- 
chymal) stem cells in vivo. Proceedings of the National Academy of Sciences 
108: 6139-6144. doi:10.1073/pnas.l016758108. 

22. Furumatsu T (2004) Smad3 Induces Chondrogenesis through the Activation of 
SOX9 via CREB-binding Protein/p300 Recruitment. Journal of Biological 
Chemistiy 280: 8343-8350. doi:10.1074/jbc.M413913200. 

23. Mar JC, Kimura Y, Schroder K, Irvine KM, Hayashizaki Y, et al. (2009) Data- 
driven normalization strategies for high-throughput quantitative RT-PCR. 
BMC Biomformatics 10: 110. doi: 10. 1 186/147 1-2105-10-1 10. 

24. Miyaki S, Nakasa T, Otsuki S, Grogan SP, Higashiyama R, et al. (2009) 
MicroRNA- 1 40 is expressed in differentiated human articular chondrocytes and 
modulates interleukin-1 responses. Arthritis Rheum 60: 2723—2730. 
doi:10.1002/art.24745. 

25. Compston J, Cooper A, Cooper C, Francis R, Kanis JA, et al. (2009) Guidelines 
for the diagnosis and management of osteoporosis in postmenopausal women 
and men from the age of 50 years in the UK. Maturitas 62: 105—108. 
doi: 1 0. 1 0 1 6/j .mattiritas. 2008. 1 1 .022. 

26. Guo L, Zhao RCH, Wu Y (201 1) The role of microRNAs in self-renewal and 
differentiation of mesenchymal stem cells. Exp Hematol 39: 608—616. 
doi:10.1016/j.cxphcm.201 1.01.011. 

27. Esau C, Kang X, Pcralta E, Hanson E, Marcusson EG, ct al. (2004) MicroRNA- 
143 Regulates .\dipocyte Differentiation. Journal of Biological Chemistry 279: 
52361-52365. doi:10.1074/jbc.C400438200. 

28. Elia L, Quintavalle M, ZhangJ, Contu R, Cossu L, et al. (2009) The knockout of 
miR- 143 and -145 alters smooth muscle cell maintenance and vascular 
homeostasis in mice: correlates with human disease. Cell Death Differ 16: 
1590-1598. doi:10.1038/cdd.2009.153. 

29. Yang B, Guo H, Zhang Y, Chen L, Ying D, et al. (2011) MicroRNA-145 
regulates chondrogenic diiferentiation of mesenchymal stem cells by targeting 
Sox9. PLoS ONE 6: e21679. doi:10.1371/joumal.pone.0021679. 

30. Lewis BP, Burge CB, Bartel DP (n.d.) TargctScanHuman. Available: http:// 
www. targe tsc an . org / cgi-bin/ targetscan/ vcrt_6 1 / targetscan. 
cgiPspecies — Human&gid — &mir_sc — miR- 1 46ac/ 1 46b-5p&mir_c — &mir_ 
nc — &mirg — &sortTypc — cs&allTxs — &incl_nc — .\11. 

31. Zhong H, Wang H-R, Yang S, ZhongJ-H, Wang T, et al. (2010) Targeting 
Smad4 links microRNA-146a to the TGF-P pathway during retinoid acid 
induction in acute promyelocytic leukemia cell Hue. Int J Hematol 92: 129—135. 
doi:10.1007/sl2185-010-0626-5. 

32. Khraiwesh B, Arif MA, Seumel GI, Ossowski S, Weigel D, et al. (2010) 
Transcriptional Control of Gene Expression by MicroRNAs. Cell 140: 111-122. 
doi:10.1016/j.cell.2009.12.023. 

33. Gao W, Yu Y, Cao H, Shen H, Li X, ct al. (2010) Deregulated expression of 
miR-21, miR- 143 and miR-181a in non small cell lung cancer is related to 
clinicopathologic characteristics or patient prognosis. Biomedicine & Pharma- 
cotherapy 64: 399-408. doi:10.1016/j.biopha.2010.01.018. 



PLOS ONE I www.plosone.org 



15 



June 2014 I Volume 9 | Issue 6 | e98063 



MiRNA-146a and Skeletal Stem Cell Differentiation 



34. Inose H, Ochi H, Kimura A, Fujita K, Xu R, ct al. (2009) A microRNA 
regulatory mechanism of osteoblast differentiation. Proceedings of the National 
Academy of Sciences 106: 20794-20799. doi:10.1073/pnas.0909311106. 

35. Klebanov L, Yalcovlev A (2007) How liigh is the level of technical noise in 
microarray data? Biol Direct 2: 9. doi:10.1 186/1745-6150-2-9. 

36. Liu Z, Lu C-L, Cui L-P, Hu Y-L, Yu Q, et al. (2011) MicroRNA-146a 
modulates T(il'-pl-induced phenotypic differentiation in human dermal 
filiroblasts by targeting SMAD4. Arch Dermatol Res 304: 195-202. 
doi:10.1007/s00403-01 1-1178-0. 

37. Jones SVV, Watkins G, Le Good N, Roberts S, Murphy CL, ct al. (2009) The 
identification of dillerentially expressed microRNA in osteoarthritic tissue that 
modulate the production of TNF-a and MMP13. Osteoarthritis and Cartilage 
17: 464-472. doi:10.1016/j.joca.2008.09.012. 

38. Liu X, Nelson A, Wang X, Kanaji N, Kim M, et al. (2009) MicroRNA- 1 46a 
modulates human bronchial epithelial cell survival in response to the cytokine- 
indueed apoptosis. Biochem Biophys Res Commun 380: 177-182. doi:10.1016/ 
j.bbrc. 2009.0 1.066. 

39. Miziino Y, rol^uzawa Y, Ninomiya Y, Yagi K, Yatsulia-Kanesalii Y, et al. (2009) 
miR-210 promotes osteoblastic differentiation tlirough inliiliition of AcvRlb. 
FEBS Lett ,i83: 2263-2268. doi:10.1016/j.iebslet.2009.06.006. 



40. Yang X, Chen L, Xu X, Li C, Huang C, et al. (2001) TGl'-beta/SmadS signals 
repress chondrocyte hypertrophic dillerentiation and are required lor maintain- 
ing articular cartilage. J Cell Biol 153: 35-46. 

41. Moustakas A, Souchelnytskyi S, Heldin CM (2001) Smad regulation in TGF- 
beta signal transduction. J Cell Sci 114: 4359-4369. 

42. Song B, Estrada KD, Lyons KM (2009) Smad signaling in skeletal development 
and regeneration. Cytokine & Growth Factor Reviews 20: 379-388. 
doi: 10.101 6/j.cytogfr.2009. 10.010. 

43. Hellingman CA, Davidson ENB, Koevoet W, Vitters EL, van den Berg WB, et 
al. (2011) Smad Signaling Determines Chondrogenic Differentiation of Bonc- 
Marrow-Derivcd Mesencliymal Stem Cells: Inhibition of Smadl/5/8P Prevents 
Terminal Differentiation and Calcification. Tissue Eng Part A 17: 1157—1167. 
doi:10.1089/trn.tra.2010.0043. 

44. Pais H, Nicolas f Iv Soond S.\I, Swinglcr TE, Clark IM, el al. (2010) .Analyzing 
mRNA expression identifies Smad3 as a microRNA- 140 target regulated only at 
protein level. RNA 16: 489-494. doi: 10.1261/rna. 17012 10. 

45. Eberhart JK, He X, Swartz ME, Yan Y-L, Song H, et al. (2008) MicroRNA 
Miml40 modulates Pdgf signaling during palatogenesis. Nat Genet 40: 290- 
298. doi:10.1038/ng.82. 



PLOS ONE I www.plosone.org 



16 



June 2014 I Volume 9 | Issue 6 | e98063 



